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ABSTRACT 

A methodology for model-independent controller 
design for controlling large angular motion of multi-body 
dynamic systems is outlined. The controlled system may 
consist of rigid and flexible components that undergo 
large rigid body motion and small elastic deformations. 
Control forces/torques are applied to drive the system, and 
at the same time suppress the vibrations due to flexibility 
of the components. The proposed controller consists of 
passive second-order systems which may be designed with 
little knowledge of the system parameters, even if the 
controlled system is non-linear. Under rather general 
assumptions, the passive design assures that the closed 
loop system has guaranteed stability properties. Unlike 
positive real controller design, stabilization can be 
accomplished without direct velocity feedback. In addition, 
the second-order passive design allows dynamic feedback 
controllers with considerable freedom to tune for desired 
system response, and to avoid actuator saturation. After 
developing the basic mathematical formulation of the 
design methodology, simulation results are presented to 
illustrate the proposed approach applied to a flexible six- 
degrcc-of- freedom manipulator. 

INTRODUCTION 

In this paper, a controller design methodology is outlined 
for non-linear dynamic systems based on simulating the 
force/torque histories that would be applied by a chosen 
virtual system of masses, springs and dashpots. In a 
companion paper by Juang and Phan*, the corresponding 
theory for linear dynamic systems is presented. Here, 
nonlinear problems are addressed by generalizing the 
standard Liapunov stability theory to allow less restrictive 
conditions on the Liapunov function. A major advantage 
of the control approach is its ability to guarantee stability 
of the controlled system with very little knowledge of the 
system being controlled. It has the advantages of positive 
real controller design methods, but in addition it allows 
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use of more realistic choices for sensors. It allows 
dynamic controllers, and it easily handles a nonlinear 
environment. A secondary advantage is the physical 
intuition available to help guide the control system 
design. Both regulator (or vibration suppression) 
problems, and end-point control (or slewing) problems 
can be treated. The approach taken here exploits the 
work-energy rate principle in designing controllers as 
suggested by Oh, Vadali and Junkins^. 

The controller design methods developed here rely on the 
guaranteed energy dissipation produced by a virtual mass- 
spring-dashpot system for their stability properties. In the 
large flexible spacecraft shape control problem, Canavin^ 
was one of the first proponents of collocated velocity 
sensors and force actuators with pure velocity feedback. 
Such a control approach has the same guaranteed stability 
properties as the proposed design method, but the 
proposed method is more general, and does not necessarily 
require velocity feedback. Hyperstability, and positive real 
controller design methods are other closely related bodies 
of knowledge, Popova-5. A positive real dynamic system 
is one which can be realized by purely passive electrical 
elements, Newcomb**. Benhabib, Iwens, and Jackson? 
were the first to suggest applying such concepts from 
network theory for shape control of large flexible 
spacecraft. Flashner^ considers applications to robotics. 
The design of linear quadratic regulator controllers that are 
constrained to be positive real is studied in Sevaston and 

Longman .9 

The proposed approach is much more general than typical 
positive real controller design. For shape control of large 
flexible spacecraft, positive real theory requires the use of 
collocated velocity sensors and force actuators (or angular 
rate sensors and torque actuators). The requirement of 
velocity sensors is very restrictive in hardware 
implementation. Here the same type of stability 
robustness properties are obtained as in positive real 
theory, but the use of position, acceleration, and velocity 
sensors is allowed, either alone or in combinations. 
Collocation of sensors and actuators is still required. The 
controller design problem is formulated in terms of a 
virtual mechanical system which supplies considerable 
intuition to guide the control system design step. 

This paper starts with a brief review of passive controller 
design for linear systems as developed in Ref. 1. The 
basic results for the linear case are then generalized to the 


1 


non-linear case. Insight into the passive control laws is 
given in the section entitled Physical Interpretation. A 
generic configuration similar to the Shuttle’s Remote 
Manipulator System (RMS) is used as an example to 
illustrate the proposed controller design approach. 

BASIC RESULTS FOR LINEAR SYSTEMS 

The virtual spring-mass-damper approach to passive 
controller design is first described for linear dynamic 
systems as in Ref. 1. Extension of the theory to general 
non-linear multiple body dynamic systems is then made. 

A general non-gyroscopic linear dynamic system can be 
represented as a system of second-order constant 
coefficient ordinary differential equations 

Mx + Dx +Kx = Bu (1) 

y = H a x + H v x + Hd x (2) 

where x is an n x 1 coordinate vector, and M, D , and K 
are the symmetric mass, damping and stiffness matrices, 
respectively. The n x p influence matrix B describes the 
actuator force distributions for the p X I control force 
vector u (or generalized force). Typically, matrix M is 
positive definite whereas D and K arc positive definite or 
positive semi-definite. In the absence of rigid-body 
motion, K is positive definite. Equation (2) is a 
measurement equation having y as the m x 1 
measurement vector, H a , // v , and Hd are the m x n 
acceleration, velocity, and displacement influence 
matrices, respectively. 

Assume that the controller to be designed has a set of 
second-order dynamic equations and measurement 
equations similar to the system equations in (1) and (2) 
above 

M c x c + D c Xc + K c Xc = B c u c (3) 

y c = Hac 'ic + live Xc + Hdc Xc (4) 

Here x c is the controller state vector of dimension n Cy and 
M c , D c , and K c are thought of as the controller mass, 
damping, and stiffness matrices, respectively. These 
matrices arc chosen to bc symmetric and positive definite 
to make the controller equations asymptotically stable. 
The n c X m influence matrix B c describes the force 
distributions for the m x 1 input force vector u c . Equation 
(4) is the controller measurement equation having y c as 
the measurement vector of length p, H ac the p x n c 
acceleration influence matrix, H vc the p x n c velocity 
influence matrix and Hdc the p x n c displacement 
influence matrix. The quantities M c , D c , K c , B Ct II ac , 
Hdo an d H vc are the design parameters for the controller. 
Let the dynamic system in (1) and (2) be connected to the 
controller system in (3) and (4) in such a way that the 
output of the controller is the input to the dynamic 


system, and tire output of the dynamic system is the input 
to the controller, i.e., 


M — yc — Hoc Xc + Hvc Xc + Hdc Xc 


Uc = y = Ha x + HvX + HdX 
The overall closed-loop system then becomes 


where 


M t xt + D t x t + K t Xt = 0 


M t = 
K t = 


M 

-BHae 
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D -BH X 

. -B c II a 

Me . 


-Belly D c . 

K 

-BHde 
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. -Be I Id 

Kc . 

L 

X c J 


(5) 

( 6 ) 

(7) 


If the system parameters, M t D, K, H a , Hd , and H v are 
known, then the controller parameters, M 0 D c , K c , H ac , 
Hdc * and H vc can be designed such that the closed-loop 
system matrices, Mt, D( y and K( are symmetric and 
positive definite. This makes the closed-loop system, Eq. 
(7), asymptotically stable. However, it is of interest to 
design controllers that are insensitive to the system 
parameters. In the following development, this is indeed 
possible provided a certain modification to the control 
equation is made, and a certain condition on the actuator 
and sensor placement is satisfied. First, note the 
similarity in structures in the closed-loop system matrices 
Mt, D( y and AT/, This means that the basic design 
procedure for the controller parameters that appear in each 
of the matrices above is the same. In particular, the case 
for designing the parameters in the stiffness matrix K t is 
illustrated here. Let the actuators be located in such a way 
that the control influence matrix of the system equation, 
B, can bc expressed by 


B T =Q b H d (8) 

then if the control influence matrix of the control 
equation, B c , is designed such that 

QbHdc = Bj (9) 

for any given matrix Hdc> then the resulting closed-loop 
stiffness matrix K t is symmetric. To derive the conditions 
that would make K ( positive definite, it is adequate to 
consider the special case where H a = H ac = H v = H vc - 0. 
Furthermore, let the input force in (5) be modified as 

u = y c - Gy = Hdc x c - GH d x (10) 

where G is a gain matrix to be determined. Also, let 

B c = K c Bc or B c = Kc l B c (11) 

and the gain matrix G be 
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G = Hdc B c 


( 12 ) 


then it can be shown that the closed-loop stiffness matrix 
in this case becomes 


K,= 


K + HdB T cK c B c Hd 

-KcBcHd 


-HdBcKc 

K c 


(13) 


T 

x a = (jc a) x a2 ••• JC^) is a generalized coordinate 

vector. The physical force or torque for each actuator is a 
function of these generalized forces through kinematic 
relations. Equation (14) is referred to as the work-energy 
rate principle prcnsented in Ref. 2. It indicates that the 
rate of change of total kinetic energy is equal to the rate of 
change of work produced by applied forces. 


which, in addition to being symmetric, is positive definite 
if the system stiffness matrix K is at least positive semi- 
definite. The design procedure for the other controller 
parameters in the closed-loop system mass and damping 
matrices are similar. 

GENERALIZATIONS TO MULTIPLE BODY 
DYNAMIC SYSTEMS 


The basic results for linear systems can be generalized to 
the non-linear case. This section outlines the overall 
methodology for extension to non-linear systems. A 
detailed development of all the possibilities involved is 
deferred to a later paper. Certain assumptions are made in 
the course of the development here which must be 
satisfied by the physical system being controlled. Various 
possible controller structures are considered for different 
types of sensor feedback, and then a Liapunov like 
approach is taken to address the stability question. First 
consider multiple rigid body systems. The case of 
multiple flexible body systems is discussed later. Also, 
mechanical systems are considered which have no external 
forces applied, such as robots in space. 

Velocity Feedback 

We start by considering controllers that use velocity 
feedback only. Since no position measurements are made, 
no attempt is made to control position. Such controllers 
are important by themselves for vibration suppression, 
and are also important here as one building block in more 
general controller designs discussed below. Let T bc the 
total kinetic energy of a mechanical system (linear or 
nonlinear) with p control actuators at p physical locations 

which arc described by p generalized coordinates Xai , i = 
1, 2, .... p. These generalized coordinates and their 
derivatives are measurable quantities such as 
displacements, velocities, accelerations or their angular 
equivalents. If the mechanical system is holonomic and 
sclcronomic (no explicit time dependence), a basic result 
of analytical mechanics relates the time derivative of the 

total kinetic energy, 7\ and the applied forces as 2 


dT 

dt 


T - 

U X n 


(14) 


Consider T as a Liapunov function 

L = L(x a ,x a ,x,x) = T(x a ,x a ,x,x) (15) 


where x denotes a vector of additional generalized 
coordinates which need not be measurable. It is assumed 
that there arc no generalized forces generated by the 
controller for these extra coordinates. The time derivative 

of the Liapunov function, L, becomes 


L 


dL 

dt 



(16) 


Let the generalized control vector u be chosen such that 


u = -Dx a 


(17) 


Then it follows that 

f—Zm. < 18 > 

dt 

which is negative semi-definite in x a ,x space when D in 
Eq. (18) is chosen as positive definite. Since the 
Liapunov function L is positive definite in the variables 

x a ,x, and its time rate of change L is only positive 
semi-definite in the same set of variables, some caution 
must be exercised with regard to a conclusion concerning 
the asymptotic stability of the overall system. If it is 
assumed that the mechanical structure of the system is 
such that x a = 0 implies 3c = 0, then it is not possible to 

have L identically zero, and yet have the kinetic energy L 
fail to vanish identically. This assumption is satisfied 

provided that the manifold of the (x 0 ,x) space on which 

the time rate of change of the Liapunov function vanishes 
identically. 


dL 

dt 




(19) 


docs not contain any arc of a trajectory, and in this case 
one can conclude that the trivial solution 


where m = (kj u 2 ••• Upf is the control vector with 
«. (/ = 1, 2, .... p) representing the generalized control 
force associated with the generalized coordinate x aii and 


x a = 0 , 3c = 0 (20) 

is asymptotically stable. It is possible, however, to 
allow the trajectories to cross the manifold, as long as 
they do not remain in it except when at the origin. The 
time rale of change of the Liapunov function considered in 
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this case is only negative semi-definite, but the fact that 
no physical trajectory arcs lie in the manifold insures that 
over any finite time interval energy is dissipated. For our 
control problem, the assumption has the physical 
implication that it is not possible to have some 
coordinates still in motion without disturbing any of the 
controlled coordinates. In practice, this condition may be 
accomplished by having a sufficient number of sensors 
collocated with the actuators at proper locations. 

Displacement Feedback 

Consider displacement feedback in order to control 
position, or alternatively to suppress vibrational motion 
as in the above section but without velocity 
measurements being available. Choose the Liapunov 
function 

L = L(x a ,x a ,x c ,x c ,x,j) (21) 

where 



where x c = (x cl x c2 ••• x cp ) is the displacement 
vector corresponding to a control (virtual) mass matrix 
M c , and M c , K c , and K are arbitrary positive definite 
matrices to be determined later. The time derivative of the 
Liapunov function, Z, becomes 

~ u *a + *c ~ *c) ^c x { x a ~~ *c ^c 2 x c 

(23) 

Let u be chosen such that 

« =~K Ct {x a ~x c ) (24) 

This is a model-independent feedback law which involves 
only the measurable generalized coordinate vector x a of 
the system and the computable displacement x c from a 
controller to be determined later. The positive definite 
matrix K c ( and the displacement vector x c are determined 

to make Z negative semi-definite. Substitution of Eq. 
(24) into Eq. (23) yields 

^ ~ “"(^a — ) Kc^X a + X c M c 'x c 

+(x fl ~~ x c | K ^ ( x a — x c ) + x c K c ^x c 
** M c x c ~ x c (r a — x c ) + x c K c ^x c 
= x T c [M e x e -K Ci {x a -x e ) + K Ci x c ] (25) 

Let the quantity in the square bracket be 


M c 'x c K C] (x a x^ + K^Xg— D c x c (26) 
or 

MX + D c x c + ( ^x c = x a (27) 

The right hand side of Eq. (27) involves only the 
measurable generalized coordinate vector x a which in turn 

determines the quantities x c , x c by computation. 
Equation (27) is a linear dynamic system represented by 
the controller mass matrix Af c , damping matrix D c and 

stiffness matrix (A^+A^) which are all positive 

definite. The generalized control vector u is generated by 
simulating the force histories that would be applied by the 
chosen virtual (controller) mass-spring-dashpot system 
described by Eq. (27). Its physical interpretation will be 
shown later. Combination of Eqs. (25) and (26) yields 

L = -x t c D c x c (28) 

which is negative semi-definite in the variables 

x a ,x a ,x ct x c ,x,x. The objective of the control is to 
obtain 

lim x a (0 = lim (0 = 0 (29) 

t~* f-»oo 

This can be accomplished provided the Liapunov function 
considered in Eq. (22) approaches zero asymptotically. 
Since L is only negative semi-definite in the variables 

x a ,x af x c ,x ct x,x, in general it is possible to have L 
vanish, yet the Liapunov function L itself does not. 
Assuming that the system is configured such that it is not 
possible to have some coordinates still in motion without 
disturbing any of the controlled coordinates, i.e., x a = 0 

implies x = 0, then in the limit L must approach zero as 
t tends to infinity. This can be seen by considering the 

case where x c = 0. Note that x c = 0 implies x a = 0. 
Therefore, the first two terms in the Liapunov function, 
Eq. (22), vanish. Furthermore, in the absence of external 

forces other than control forces, assume that x c = 0 
implies x c =0, which from Eq. (27) yields x fl =0 

provided K Ci is invertible. This makes the last two 
terms in the Liapunov function vanish as well. 

Computation of the quantity x c requires measurement of 
the vector x a as shown in Eq. (27). Therefore, when 
velocity sensors arc not available, so that the control law 
cannot depend on measuring x a as in the previous 

section, it is still possible to make L negative semi- 
definite using the feedback law given in Eq. (24), and 
accomplish vibration suppression The positive definite 

matrices M c% D c , and K ^ are arbitrary and can be 
chosen to meet performance requirements. 
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Displacement and Acceleration Feedback 

Feedback laws using acceleration measurements can be 
derived similarly. This is an important case for vibration 
suppression in structures since acceleration measurements 
are much more easily made than velocity measurements. 
Consider a Liapunov function of the form 

L = T + ±(i. + x c f M c {x a + i c )+LX* f 

+ l(x a + x e ) T K Ci (x a + x c ) (30) 

The time derivative of L then becomes 

= u T x a + (x a + x c ) T M c [x a + x c ) + xj K c x c 
at 

+ {x a + x c ) T K Ci (x a + x c ) ( ' 1} 

Let the control input be chosen such that 

u = -M c {x a + x c )-K Ci {x a + x c ) (32) 

then the time rate of change of the Liapunov function 
becomes 


measurements is sufficient to make L negative semi- 
definite. Again, assuming that x a =0 implies xsO i it 
can be concluded that L approaches zero as t tends to 
infinity. The line of argument follows closely the 
previous development in the case of displacement feedback 
and hence is omitted. 

Displacement, Velocity, and Acceleration 
Feedback 

The general case of displacement, velocity, and 
acceleration feedback can be easily derived by considering 
a Liapunov function of the form 

L = T + L X +L 2 (37) 

where 

L 1 ~ ^ X cd^ cd X cd 

1 t It 

+ ” (^ fl ~ X c 4 ) { x a — X cd ) ^ 2^ x cd^cd l x cd 

1 T 

L 2 = { x a + x ca ) ^ ca (*a + X ca ) 

1 1 T 

+ ^ x ca^ca [ X ca + ^ ( X a + X ca ) ( X a + X ca ) 


^ = -(* a + X c f M c x a - {x a + f K c x a 
at 


+(i a + x c ) M c + x c ) + x c K c x c 
= x^M c [x a + x c ) + x c K c ) x c + x c K c ^ (x a + x c ) 

K 

- { X a + f K cSa + £ (*• + ** ) 

1 = 1 

= x c Cj x c + M c (x a + x c ) + K c ^ ( x a + x c ) j 


(33) 


If the quantity in the square bracket in the above equation 
is set to be 


Kc*c + M c {x tt + x c ) + K C2 {x a + x c ) = -D c x c (34) 
or 

M c x c + D c x c + (K Cx + K h )x e = -M c x a - K H x a (35) 

where D c is an arbitrary positive definite matrix, then Eq. 
(33) becomes 

2—gDA <36) 

dl 


which is negative semi-definite in the variables 

x tt ,x a ,x c ,x c ,x,t. Note that even in the absence of 
velocity measurements, the feedback control law given in 
Eq. (32) which is based on acceleration and displacement 


The control law is then of the form 

U = —Dx a — — %cd ) — M ca {x a + X ca ) 

~K cai (x a + x ca ) (38) 

where the vectors X cd and X ca are determined from 

Media + D aia + { K a, + K cd 2 ) x a = ~ K cd > x a 
M ca i ca + D ca x ca + jx ca — ~M ca x a — K ca x x a 

The coefficient matrices as well as K cJt and K c<h in Eq. 

(39) are all required to be positive definite. The time rate 
of change of the Liapunov function then becomes 

L = -igDx a — x c jD cc ix c 4 ~ x ca D ca x ca (40) 

which is negative semi-definite. The proof of stability 
parallels the previous case. 

Generalization to Multiple Flexible Body 
Systems with Viscous Friction 

Including flexibility into the bodies of the system, 
introduces a potential energy of deformation. In many 
cases such as robots with controllers at each joint, this 
potential energy is special in the sense that it does not 
alter the equilibrium position of the system including the 
virtual controller. Since the potential energy function does 
not change the equilibrium position, or introduce new 
equilibria, there is no need to solve for the equilibrium 
before studying stability. To include such a potential 
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energy in the analysis, one simply notes that the time rate 
of change of the kinetic energy in Eq. (14) becomes the 
time rate of change of the total energy *• so that 

±(T+V)=-2F d + u T x a (41) 

at 

where F d is the quadratic Rayleigh's dissipation function 
included here to handle any viscous damping in the 
system. All the theory developed above applies again. 
One simply replaces the kinetic energy T by the total 
mechanical energy E = T + Vof the system in the related 
expressions. 

PHYSICAL INTERPRETATION 

Insight into the control laws developed can be gained 
through the following physical interpretation in terms of 
mechanical analogs. It is these analogs that motivated the 
mathematical formulations above. For simplicity, 
consider only the single input and single output case. 

The matrices M c , D c , K c reduce to scalar quantities and 

are denoted by the lower case letters m c , d c> k c . 

Velocity Feedback 

The controller is given by Eq. (17). The control force 
applied is equivalent to the force transmitted to the system 
at the actuator location by a dashpot with a damping 
coefficient d. This is shown in Fig. 1. 



Xa 

Figure 1: Velocity Feedback 

If X a is an angular displacement, then the control law 
corresponds to an angular dashpot 


Displacement Feedback 

The control law is given in Eq. (24). Examination of 
Eqs. (24) and (27) shows that the control force is 
equivalent to that provided by a mass-spring-dashpot 
system at the actuator location. This is shown in Fig. 2. 


trie dc 



Xa Xc 


Figure 2: Displacement Feedback 

In this case, the quantity x a is the measured inertial 
position of the system at the actuator location, whereas 
the quantity x c is defined as the inertial position of the 

virtual mass m c . The springs are considered as ideal 
springs, and their unstrelched lengths can be considered to 
be zero. The corresponding block diagram for this case is 
given in Fig. 3. 

Displacement and Acceleration Feedback 

The control law is given in Eqs. (32) and (35). 
Examination of the controller equations verifies that the 
control force is equivalent to a mass-spring-dashpot 
system at the actuator location, with one important 
distinguishing feature. The position of the controller 

(virtual) mass x c shown in Fig. 4 is now defined relative 
to the measured inertial position x a of the system at the 
actuator location whereas in the displacement feedback 



Figure 3: Displacement Feedback Controller 
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case above, it is defined as an inertial position. The free 
lengths of the springs are taken as zero. The 
corresponding block diagram for this case is given in Fig. 


k c , 


m c 


^ | — — 


K 

V 


kci 


Xa 


Figure 4: Displacement and Acceleration 
Feedback 

Displacement , Velocity, and Acceleration 
Feedback 

The control law is given in Eq. (38). The equivalent 
mass-spring-dashpot diagram for this case is given in Fig. 
6 and the block diagram for this controller is given in 
Fig. 7. 


Xa 


-wmM 
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I — VvVWvV — 
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Figure 6: Displacement, Velocity, and 
Acceleration Feedback 



Figure 5: Displacement and Acceleration Feedback Controller 



Figure 7: Displacement, Velocity, and Acceleration Feedback Controller 
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REMARKS 


Remark 1: Certain simplifications of the controllers 
developed here are worth noting. In both the displacement 
feedback, and displacement and acceleration feedback cases, 
the inclusion of the spring characterized by its coefficient 

is necessary for position control, but not vital for 

vibration suppression. This special feature allows one to 
design robust vibration suppression controllers with 

accelerometer feedback. In particular, setting k^ = 0 in 

the displacement and acceleration feedback law, Eq.(32), 
results in a feedback law that is based on acceleration 
measurements alone. The corresponding diagram for this 
case is given in the Fig. 8. 


^ j — — 


Figure 8: Acceleration Feedback 

Remark 2: In the absence of direct velocity 

measurements, it is possible to mimic direct velocity 
feedback by displacement measurements alone. Consider 
the case of displacement feedback as shown in Fig. 2. By 

letting k H ~ 0, and choosing k c% to be large, and m c to 

be small, the behavior of the resulting mass-spring- 
dashpot system approximates that of a single dashpot. 
This is equivalent to putting a controller zero near the 
origin but not at the origin. 

Remark 3: It is also possible to approximate direct 
velocity feedback from displacement and acceleration 

measurements. Refer to Fig. 4. If k q , m c are chosen to 

be small, and to be large, then the behavior of the 

resultant mass-spring-dashpot system approximates that 
of a single dashpot. In the absence of displacement 

measurement ( k = 0), the same effect can be achieved by 
choosing k Ci to be small, but m c to be large. 

Remark 4: To implement the controller in Fig. 2 or 4, 
it is necessary to select the correct initial conditions for 
the controller coordinate x c to be used. There is a subtle 
difference between the vibration suppression problem and 
the end-point control problem by this approach, even 
though they may be viewed as the same conceptually. 
For vibration suppression, the controller consists of 
second-order dynamic systems whose equilibrium position 
of the controller mass is consistent with the desired 


equilibrium configuration of the structure. Assuming that 
the structure is initially at rest, the initial condition for 

the controller coordinate is then x c - x c = 0. For end- 
point positioning control, however, the initial conditions 
for the controller coordinate must be chosen such that in 
the desired final configuration of the manipulator, the 
virtual spring and dashpot are in their equilibrium and 
unstretched position. In the end-point control problem, 
the induced vibrations are also suppressed by this 
controller design. 

Remark 5: When the controllers are interpreted in 
physical terms, recall that for the case of displacement 
feedback, the position vector of the controller (virtual) 

mass x c in the Liapunov function, Eq. (22), is defined 
with respect to an inertial reference. However, for 
displacement and acceleration feedback, Eq. (30), it is 
defined relative to the actuator location. When the control 
mass location is expressed as an inertial position, the 
Liapunov function in the second case becomes identical to 
that in the first case. The difference is in the formulation 
of the control law. By expressing it in terms of the 
relative position of the controller mass, the scheme can be 
implemented by acceleration measurements without the 
need for velocity measurements. 

NUMERICAL EXAMPLES 

The robot system, shown in Fig. 9, consists of six one 
degrec-of-freedom joints, namely, shoulder pitch, yaw, 
elbow pitch, wrist yaw, pitch and roll. This robot model 
has the same kinematic relationships and mass properties 
as the Shuttle's Remote Manipulator System (RMS). 
The shoulder and elbow joints are connected by a 6,4 
meter long lightweight carbon composite boom. This is 
designated the upper arm boom. The lower arm boom, 
connecting the elbow joint to the wrist joints, is 
approximately 7 meters long. The shoulder and elbow 
joints provide three translational degrees of freedom, and 
the wrist joints provide three additional rotational degrees 
of freedom of the end effector. Mass properties of each of 
the components of the robot system are given in Table 1. 

A simple maneuver to move the end effector from its 
current position to a new position is simulated. Initially, 
the robot system is in the position shown in Fig. 9, and 
the position vector of the end effector is (15.3162, 

0. 3048, 0.0) given in terms of its x, y, z components. 
The desired new position of the end effector relative to the 
robot base is (10.3162, 5.3048, 5.0). Assume that all the 
wrist joints are locked during the maneuvering, and 
control torques arc only applied to each movable joint, 

1. e., the shoulder yaw, pitch, and elbow pitch joints. 
Since the system is driven by torques at the joints, the 
required change in each joint angle is pre-calculated to be 
-0.45131 rad,-0.3111 rad, and 1.2106 rad for the shoulder 
yaw, pitch, and elbow pitch joints, respectively. In this 
case, a second-order virtual rotational spring-mass damper 
system is located at each joint to be controlled. The 
controller system then has virtual moments of inertia 
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instead of virtual masses* torsional springs and dashpols 
instead of linear springs and dashpots. The initial 
conditions of the control system arc chosen such that in 
the final desired configuration of the manipulator, the 
virtual springs assume their free-length forms. 

In the simulations, the wrist links are modeled as rigid 
bodies since they are fairly short relative to the lengths of 
both arms. The upper and lower arms are modeled as 
flexible bodies. The geometric and material properties of 
the flexible arms are given in Table 2. Only the first three 
cantilever modes are used to represent the flexibility of the 
links. Three different controllers, distinguished by the 
measured signals, are used to illustrate the use of the 
controller design methods in this paper. 

Case 1: Displacement and velocity feedback 

Case I uses angular displacements and velocities of joints 
for feedback. The controllers for each joint are rotational 
versions of the diagram in Fig. 6 with the bottom 
connection associated with acceleration removed. The 
parameters of the controller are given in Table 3. No 
attempt has been made to optimize the parameters chosen. 
Figures 10-12 show time histories of the end effector 
locations in both flexible and rigid-body simulations. The 
results show that the vibration is actively suppressed. 
Figure 13 shows the time history of the dominant first 
mode of the lower arm during this maneuver. 

Case 2: Displacement feedback 

One advantage of the virtual system controller design 
approach developed here is that velocity feedback is not 
required for vibration control or stability robustness. The 
stability robustness however, need not guarantee good 
damping rates. To auempt to produce responses like those 
in Case 1, but this time using displacement measurement 
feedback alone, the order of the controller is increased to 
emulate velocity feedback. This is done by replacing the 
velocity feedback portion of Case 1 with a rotational 
virtual system corresponding to Fig. 2, and using a large 

k C] and a small m c . The control structure can again be 

viewed by comparing to Fig. 6. The bottom acceleration 
feedback connection is removed as before, and in place of 
the middle velocity connection, we substitute a repeat of 
the top connection, but with a different choice of 

parameters. Theoretically, when both the k c ^ and m c of 

the middle connection go to their limits, the control 
system in Case 2 is equivalent to the control system in 
Case 1. Numerically, the values of these constants are 
constrained by the introduction of high frequencies and 
real-time numerical integration difficulties. 

In the simulation, the value of k^ of the second virtual 

system is taken to be 100 times suffer than the k Cx of the 
first virtual system. Figure 14 shows the simulation 
results with two different values of m c for the second 


virtual system, and shows that this approach will require 
rather extreme mass and stiffness values for good damping 
performance. 

Case 3: Displacement and displacement-acceleration 
feedback 

Another way to emulate the velocity signal is by using 
displacement-acceleration feedback. By replacing the 
velocity feedback in Case 1 with the displacement- 

acceleration feedback of Fig. 4 and using large k^ and 

small m c in this virtual system, the controller works 
similarly to the controller used in Case 1. It can be 
viewed as a rotational version of Fig. 6 with the center 

velocity connection removed. Theoretically, when the k ^ 

and m c of the acceleration connection go to their limits, 
the control system in Case 3 is equivalent to the control 
system in Case 1. The controller in Case 3 is numerically 
better conditioned than that in Case 2, since the velocity 
signal is obtained in the limit through integration of 
acceleration rather than differentiation of displacement 

In the simulation, the value of k^ of the second virtual 

system is taken to be 100 times suffer than the k Ci of the 

first virtual system. The results obtained are identical to 
within plotting accuracy, to the results obtained in Case 
1, and are therefore not shown here. Since acceleration 
feedback is often more realistic than velocity feedback, 
this example shows that the control design of Case 3 can 
be important in practical applications. 

So far, the controller parameters have not been optimized. 
To get some understanding of how the parameters of the 
virtual system affect performance, several cases are 

simulated with adjusted m c , k c , and d c of the first 

virtual control system in Case 3. Each of the parameters 
is changed to ten times larger and ten times smaller than 
its original value. Figure 15 shows the angular 
displacement of the elbow pitch joint, for different values 

of k Ci . Faster response can be achieved by using a larger 
spring constant k Cj , as the results revealed. Figure 16 

gives results obtained when varying the control mass m c t 

In this particular example, the smaller the controller 
mass, the better the response. Figure 17 presents the 

results of varying the control damping constant d c . Note 
that the use of a large damping constant may result in 
slow system response. 

CONCLUSIONS 

In this paper, a general control design methodology has 
been outlined for large angle position control and 
vibration suppression in multiple flexible-body dynamic 
systems. The method guarantees stability of the 
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controlled system, whether linear or non-linear. The 
approach is model-independent in the sense that 
knowledge of the system dynamics is not required in the 
design process. Hence, it is robust with respect to 
parameter variations. Unlike positive real controller 
design which allows velocity measurements only, the 
current development can use velocity, position, and 
acceleration measurements or any combination, and still 
guarantees stability of the closed loop system. In addition 
it allows use of dynamic feedback controllers with 
additional opportunity for tuning to obtain good system 
performance. One special case of the methodology 
demonstrates how one can obtain the desirable properties 
of velocity feedback using acceleration feedback, which is 
often easier to instrument in practice. The design has 
intuitive appeal in terms of its physical interpretation that 
aids the control design process. 

REFERENCES 

* JUANG, J.-N. and PHAN, M., "Robust Controller 
Designs for Second-Order Dynamic Systems: A Virtual 
Passive Approach," NASA Technical Memorandum TM - 
102666 , May 1990. 

2 OH, H-S., VADALI. S.R., and JUNKINS, J.L., "Use 
of the Work-Energy Rate Principle for Designing 
Feedback Control Laws," Journal of Guidance, Control , 
and Dynamics, to appear. 


3 CANAVIN, J.R., "The Control of Spacecraft Vibration 
Using Multivariable Output Feedback," Proceedings of 
the AIAAIAAS Astrodynamics Conference, Paper No. 78- 
1419, Palo Alto, California, August 1978. 

4 POPOV, V.M., Hyperstability of Automatic Control 
Systems, Springer- Verlag, New York, 1973. 

5 POPOV, V.M., Lectures on Stability Theory, 
University of Southern California, Electrical Engineering 
Technical Report 357, 1969. 

6 NEWCOMB, R.W., Linear Multiport Synthesis, 
McGraw-Hill, New York, 1966. 

7 BENHABIB R.J., IWENS, R.P., and JACKSON, R.L., 
"Stability of Large Space Structure Control Systems 
Using Positivity Concepts," Journal of Guidance, 
Control, and Dynamics, Vol. 4, No. 5, Sept.-OcL 1981. 

* FLASHNER, H„ "An Approach for Control of Robot 
Manipulators," Proceedings of the AIAA Guidance, 
Navigation, and Control Conference, Snowmass, 
Colorado, August 1985. 

9 SEVASTON, G. and LONGMAN, R.W., "Optimal 
Positive Real Controllers for Large Space Structures," 
Proceedings of the Third VP1&SUIAIAA Symposium on 
Dynamics and Control of Large Flexible Spacecraft, L. 
Meirovitch, Editor, Blacksburg, Virginia, June 1981. 


10 
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Tnble 2 Geometries of the Robot System and Material 
Properties of Flexible Links 
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Table 3 Parameters of Virtual Systems for the Robot Maneuvering 
Control with Velocity and Displacement Feedback 
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Figure 9: A Six-Degree-of-Freedom Robot 



Figure 10: Time Histories of the x-Position 
of the End Effector, Case 1 
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MODAL DISPLACEMENT 



Figure 12: Time Histories of tire z-Position 
of the End Effector, Cnse 1 



Figure 13: Time Histories of the First Mode 
of the Lower Arm during Maneuver 



Figure 14: Time Histories of the x-Posidon 
of die End Effector, Case 2 
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Figure 15: Time Histories of the Displacement of the Elbow Pitch 
for Different Values of Control Stiffness Constant 



Figure 16: Time Histories of the Displacement of the Shoulder Pitch 
for Different Values of Control Mass Constant 



Figure 17: Time Histories of the Displacement of the Shoulder Yaw 
for Different Values of Control Damping Constant 
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